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Abstract 

The production of prompt D+ mesons was measured for the first time in collisions of heavy nuclei 
with the ALICE detector at the LHC. The analysis was performed on a data sample of Pb-Pb 
collisions at a centre-of-mass energy per nucleon pair, of 2.76 TeV in two different centrality 

classes, namely 0-10% and 20-50%. D+ mesons and their antiparticles were reconstructed at 
mid-rapidity from their hadronic decay channel D+ —> with 0 —i>K^K+, in the transverse 

momentum intervals A < pj < \2 GeV/c and 6 < pi < 12 GeV/c for the 0-10% and 20-50% 
centrality classes, respectively. The nuclear modification factor /?aa was computed by comparing 
the pT-differential production yields in Pb-Pb collisions to those in proton-proton (pp) collisions at 
the same energy. This pp reference was obtained using the cross section measured at ^/s = 1 TeV and 
scaled to i/s = 2.76 TeV. The /?aa of D+ mesons was compared to that of non-strange D mesons 
in the 10% most central Pb-Pb collisions. At high pj (8 < px < 12 GeV /c) a suppression of the 
-meson yield by a factor of about three, compatible within uncertainties with that of non-strange 
D mesons, is observed. At lower px (4 < px < 8 GeV /c) the values of the D+-meson 7 ?aa are larger 
than those of non-strange D mesons, although compatible within uncertainties. The production ratios 
D+ /DO and D+ /D+ were also measured in Pb-Pb collisions and compared to their values in proton- 
proton collisions. 
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1 Introduction 

Calculations of Quantum Chromodynamics (QCD) on the lattice predict that strongly-interacting matter 
at temperatures exceeding the pseudo-critical value of about ~ 145 — 165 MeV and vanishing baryon 
density behaves as a deconfined Plasma of Quarks and Gluons (QGP) mill. In this state, paitons are the 
relevant degrees of freedom and chiral symmetry is predicted to be restored. The conditions to create a 
QGP are expected to be attained in collisions of heavy nuclei at high energies. This deconfined state of 
matter exists for a short time (few fm/c), during which the medium created in the collision expands and 
cools down until its temperature drops below the pseudo-critical value and the process of hadronisation 
takes place. 

Heavy quarks (charm and beauty) are sensitive probes to investigate the properties of the medium formed 
in heavy-ion collisions. They are produced in quark-antiquark pairs predominantly at the initial stage 
of the collision in hard-scattering processes characterized by timescales shorter than the QGP formation 
time ina. The heavy quarks propagate through the expanding hot and dense medium, thus experiencing 
the effects of the medium over its entire evolution. While traversing the medium, they interact with its 
constituents via both inelastic and elastic QCD processes, exchanging energy and momentum with the 
expanding medium iHm. For heavy quarks at intermediate and high momentum, these interactions lead 
to energy loss due to medium-induced gluon radiation and collisional processes. 

Evidence for heavy-quark in-medium energy loss is provided by the observation of a substantial modifi¬ 
cation of the transverse momentum (px) distributions of heavy-flavour decay leptons iT/l- fTOl . D mesons 
II11I12I1 and non-prompf J/t/r |[T3l in Au-Au and Pb-Pb collisions af RHIC and EHC energies as compared 
fo profon-profon (pp) collisions. This modification is usually quantified by fhe nuclear modification fac- 
for Raa> defined as fhe ratio befween fhe yield measured in nucleus-nucleus collisions and fhe cross sec- 
fion in pp inferacfions scaled by fhe average nuclear overlap function. In absence of nuclear effecls, Raa 
is expecfed fo be unify. Parfon in-medium energy loss causes a suppression of hadron yields, /?aa < 1> 
af intermediate and high fransverse momenfum (px > 3 GeV/c). In cenfral nucleus-nucleus collisions 
af RHIC and EHC energies, Raa values significanfly lower fhan unify were observed for heavy-flavour 
hadrons wifh px values larger fhan 3—4 GeV/c. In fhis px range, fhe D-meson yields measured in p-Pb 
collisions af = 5.02 TeV are consisfenf wifh binary-scaled pp cross sections l[T4ll . providing clear 
evidence fhaf fhe suppression observed in Pb-Pb collisions is nof due fo cold nuclear mailer effecls and 
is induced by a slrong coupling of fhe charm quarks wifh fhe hoi and dense medium. 

In case of subslanfial inferacfions wifh fhe medium, heavy quarks lose a significanf amounf of energy 
while Iraversing fhe fireball and may participate in fhe collective expansion of fhe system and possibly 
reach Ihermal equilibrium wifh fhe medium consfiluenfs. In fhis respecl, fhe measuremenl of a positive 
elliptic flow V 2 of D mesons af EHC energies ifTSlflfill and of heavy-flavour decay eleclrons af RHIC 
energies HHUllTl provides an indicalion fhaf fhe inferacfions wifh fhe medium consfiluenfs fransfer fo 
charm quarks informalion on fhe azimufhal anisofropy of fhe sysfem. 

If is also predicted fhaf a significanf fraction of low- and inlermediale-momenlum heavy quarks could 
hadronise via recombination wifh ofher quarks from fhe medium |[T8l420l . An imporlanf role of hadro- 
nisalion via (re)combinalion, eifher during fhe deconfined phase 11211 or af fhe phase boundary l[22ll . is 
indeed supporfed by fhe resulls of J/t/r nuclear modificalion facfor and elliptic flow af low px l[23ll25]l . 
Hadronisalion via recombinalion allows in some models, e.g. |[26] - l28l . a better descriplion of heavy- 
flavour production measuremenls af RHIC and EHC energies, in particular fhe /?aa of D** mesons af low 
Px measured in Au-Au collisions af = 200 GeV l[T2ll and fhe positive and sizable D-meson V 2 in 
Pb-Pb collisions af = 2.76 TeV ifTSll . 

The measuremenl of D+-meson production in Pb-Pb collisions can provide crucial addilional informa¬ 
lion for underslanding fhe inferacfions of charm quarks wifh fhe slrongly-inferacling medium formed 
in heavy-ion collisions af high energies. In parlicular, fhe D+-meson yield is sensitive fo sfrangeness 
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production and to the hadronisation mechanism of charm quarks. 

An enhancement of strange particle production in heavy-ion collisions as compared to pp interactions 
was long suggested as a possible signal of QGP formation II291I30II . Strange quarks are expected to be 
abundant in a deconfined medium due to the short time needed to reach equilibrium values among the 
parton species and to the lower energy threshold for ss production. A pattern of strangeness enhancement 
increasing with the hadron strangeness content when going from pp (p-A) to heavy-ion collisions 
was observed at the SPS EIHll, at RHIC 133 and at the LHC In the frame of the statistical 
hadronisation models, strange particle production in heavy-ion collisions follows the expectation for a 
grand-canonical ensemble. In contrast, for pp collisions canonical suppression effects are found to be 
important, reducing the phase space available for strange particles II371I38II . In this context, the increase 
in strange particle yields in heavy-ion collisions compared to pp interactions is viewed as due primarily 
to the lifting of the canonical suppression. 

This strangeness enhancement effect could also affect the production of charmed hadrons if the dominant 
mechanism for D-meson formation at low and intermediate momenta is in-medium hadronisation of 
charm quarks via recombination with light quarks. Under these conditions, the relative yield of D+ 
mesons with respect to non-strange charmed mesons at low pj is predicted to be enhanced in nucleus- 
nucleus collisions as compared to pp interactions 11391 - 1411 . The comparison of the pT-differential 
production yields of non-strange D mesons and of D+ mesons in Pb-Pb and pp collisions is therefore 
sensitive to the role of recombination in charm-quark hadronisation. 

A consequence of the possibly enhanced production of D+ mesons in heavy-ion collisions would be a 
slight reduction of the fraction of charm quarks hadronising into non-strange meson species. Therefore, 
the measurement of the D+-meson production is also relevant for the interpretation of the comparison 
of the nuclear modification factors of non-strange D mesons and light-flavour hadrons (pions) 1111 [1421 . 
which is predicted to be sensitive to the quark-mass and colour-charge dependence of paiton in-medium 
energy loss ||3|43j|44l. Furthermore, due to this possible modification of the relative abundances of D- 
meson species, measuring the D+ yield at low pj is needed also to determine the total charm production 
cross section in Pb-Pb collisions. 

The px-differential inclusive production cross section of prompQ D+ mesons (average of particles and 
antiparticles) was measured in pp collisions at ^/s = 7 TeV with the ALICE detector and it was found 
to be described within uncertainties by perturbative QCD (pQCD) calculations 1451 . The D+ nuclear 
modification factor was measured in p-Pb collisions at = 5.02 TeV and found to be consistent 
with unity 11411 . In this paper, we report on the measurement of prompt D+-meson production and 
nuclear modification factor in Pb-Pb collisions at = 2.76 TeV. D+ mesons (and their antiparticles) 
were reconstructed at mid-rapidity, |y| < 0.5, through their hadronic decay channel D+ —)■ ^7r+ with a 
subsequent decay 0 —)• K^K+. The production yield was measured in two classes of collision centrality, 
central (0-10%) and semi-central (20-50%), and compared to a binary-scaled pp reference obtained by 
scaling the cross section measured at ^/s = 1 TeV to the Pb-Pb centre-of-mass energy via a pQCD-driven 
approach. The experimental apparatus and the data sample of Pb-Pb collisions used for this analysis are 
briefly presented in Section [2l In Section [3l the D+ meson reconstruction strategy, the selection criteria 
and the raw yield extraction from the KKtt invariant mass distributions are discussed. The corrections 
applied to obtain the pj-differential production yields of D+ mesons, including the subtraction of the 
non-prompt contribution from beauty-hadron decays, are described in Section ID The various sources of 
systematic uncertainty are discussed in detail in Section[5l The results on the D+-meson production yield 
and nuclear modification factor are presented in Section |6] together with the comparison to non-strange 
D-meson /?aa and to model calculations. The D+/D® and D+/D+ yield ratios in three pj intervals for 

*In this paper, ’prompt’ indicates D mesons produced at the interaction point, either directly in the hadronisation of the 
charm quark or in strong decays of excited charm resonances. The contribution from weak decays of beauty hadrons, which 
gives rise to feed-down D mesons displaced from the interaction vertex, was subtracted. 
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the 10% most central Pb-Pb collisions are compared to those in pp collisions. 

2 Apparatus and data sample 

The ALICE detector and its performance are described in detail in Refs, ll^ and HTll . respectively. 
The apparatus consists of a central barrel covering the pseudorapidity region |t] | < 0.9, a forward muon 
spectrometer (—4.0 < rj < —2.5) and a set of detectors for triggering and event centrality determination. 
The detectors of the central barrel are located inside a 0.5 T magnetic field parallel to the LHC beam 
direction, that corresponds to the z-axis in the ALICE reference frame. The information provided by 
the following detectors was utilised to perform the analysis presented in this paper: the Inner Tracking 
System (ITS), the Time Projection Chamber (TPC) and the Time Of Plight (TOP) detector were used 
to reconstruct and identify charged particles at mid-rapidity, while the VO scintillator detector provided 
the information for triggering, centrality determination and event selection. The neutron Zero Degree 
Calorimeters (ZDC) were also used, together with the VO detector, for the event selection. 

The trajectories of the D-meson decay particles are reconstructed from their hits in the ITS and TPC 
detectors. Particle identification is performed utilising the information from the TPC and TOP detectors. 
The ITS consists of six cylindrical layers of silicon detectors covering the pseudorapidity interval 
|t]| < 0.9. The two innermost layers, located at 3.9 and 7.6 cm from the beam line, are composed of 
Silicon Pixel Detectors (SPD). The two intermediate layers are equipped with Silicon Drift Detectors 
(SDD) and the two outermost layers, with a maximum radius of 43.0 cm, are composed of double-sided 
Silicon Strip Detectors (SSD). The high spatial resolution of the ITS detectors, together with the low 
material budget (~ 7.7% of a radiation length at T] = 0) and the small distance from the interaction point, 
provides a resolution on the track impact parameter (i.e. the distance of closest approach of the track to the 
primary vertex) better than 65 /im for transverse momenta pt > 1 GeV/c in Pb-Pb collisions BTl . The 
TPC, covering the pseudorapidity interval |t]| < 0.9, provides track reconstruction with up to 159 points 
along the trajectory of a charged particle and allows its identification via the measurement of specific 
energy loss dL/dx. Particle identification is complemented with the particle time-of-llight measured with 
the TOP detector, which is composed of Multi-gap Resistive Plate Chambers and is positioned at 370-399 
cm from the beam axis, covering the full azimuth and the pseudorapidity interval |t]| < 0.9. The TPC 
and TOP information provides pion/kaon separation at better than 3 a level for tracks with momentum 
up to 2.5 GeV/c HTl . 

The analysis was performed on a sample of Pb-Pb collisions at centre-of-mass energy per nucleon pair, 
of 2.76 TeV collected in 2011. The events were recorded with an interaction trigger that required 
coincident signals in both scintillator arrays of the VO detector, covering the pseudorapidity ranges 
—3.7<T]<—1.7 and 2.8 < q < 5.1, respectively. An online selection based on the VO signal amplitude 
was used to record samples of central and semi-central collisions through two separate trigger classes. 
Events were further selected offline to remove background from beam-gas interactions on the basis of 
the timing information provided by the VO and the neutron ZDC detectors (two hadronic calorimeters 
located at z = 114 m on both sides of the interaction point covering the interval |t] | > 8.7). Only events 
with an interaction vertex reconstructed from ITS-i-TPC tracks with |zvertex| <10 cm were considered in 
the analysis. 

Collisions were classified in centrality classes based on the sum of the signal amplitudes in the two 
VO scintillator arrays. Each class is defined in terms of percentiles of the hadronic Pb-Pb cross 
section, as determined from a fit to the VO signal amplitude distribution based on the Glauber-model 
description of the geometry of the nuclear collision ll48]450ll and a two-component model for particle 
production lISll . The analysis was performed in two centrality classes: 0-10% and 20-50%. In total, 
16.5 xlO^ events, corresponding to an integrated luminosity Lint = (21.5 ±0.7) /ub \ were analysed in 
the 0-10% centrality class, and 13.5x10^ events. Lint = (5.9±0.2) jUb in the 20-50% class. The 


4 



Ds production in Pb-Pb collisions 


ALICE Collaboration 


average values of the nuclear overlap function Taa (defined as the convolution of the nuclear density 
profiles of the colliding ions |[50ll and proportional to the number Acoii of binary nucleon-nucleon 
collisions occurring in the Pb-Pb collision) are reported in Table[T]for the 0-10% and 20-50% centrality 
classes, together with their systematic uncertainty estimated as described in lISTl . 


Cenfralify class 

(Taa) (mb ^) 

Aevt 

^int (/^b 

0-10% 

23.44 zb 0.76 

16.4x10^ 

21.3zb0.7 

20-50% 

5.46 ± 0.20 

13.5x10^ 

5.8±0.2 


Table 1: Average value of the nuclear overlap function, (7aa), for the considered centrality classes, expressed as 
percentiles of the hadronic Pb-Pb cross section. The values were obtained with a Monte Carlo implementation of 
the Glauber model assuming an inelastic nucleon-nucleon cross section of 64 mb ED. The number of analysed 
events and the corresponding integrated luminosity in each centrality class are also shown. The uncertainty on 
the integrated luminosity derives from the uncertainty of the hadronic Pb-Pb cross section from the Glauber 
model lISTI . 


3 D+ meson reconstruction and selection 

D+ mesons and their antiparticles were reconstructed in the decay channel D+ —)• 071+ —)■ K^K+7r+ (and 
its charge conjugate), whose branching ratio (BR) is (2.24 zb 0.10)% ESI. Other D+ decay channels can 
give rise to the same K^K+7r+ final sfafe, such as D+ —K***K+ and D+ —)■ fo(980)7r+, wifh BR of (2.58 
zb 0.11)% and (1.14 zb 0.31)%, respecfively Il5^ . However, as explained in Ref. Il45]l . fhe applied cufs for 
fhe selection of fhe D+ signal candidafes sfrongly reduce confribufions from fhese channels, and Iherefore 
fhe measured yield is dominated by fhe D+ —)• 071+ —)■ K^K+7r+ decays. The decay channel fhrough fhe 
0 resonance was chosen because fhe narrower widfh of fhe 0 invarianf-mass peak with respect to fo (980) 
and K provides the best discrimination between signal and background. 

The analysis strategy for the extraction of the signal out of a large combinatorial background is based 
on the reconstruction of decay topologies with a secondary vertex significantly displaced from the 
interaction point. The secondary vertex position and its covariance matrix were determined from the 
decay tracks by using the same analytic minimization method as for the computation of the primary 
vertex 1531. The resolution on the position of the D+ decay vertex was estimated with Monte Carlo 
simulations and it was found to be about 100 /rm. D+ mesons have a mean proper decay length 
cT = 150 zb 2 /rm l52l . which makes it possible to resolve their decay vertices from the primary vertex. 
With the current data sample, the signal of D+ mesons could be extracted in three pj intervals (4-6, 6-8 
and 8-12 GeV/c) in the 0-10% centrality class and in two pj intervals (6-8 and 8-12 GeV/c) in the 
20-50% centrality class. 

D+ candidates were defined from friplefs of fracks wifh fhe proper charge sign combination. Tracks 
were selected requiring |t]| < 0.8 and px > 0.6 (0.4) GeV/c in fhe 0-10% (20-50%) cenfralify class. In 
addition, fracks were also required fo have af leasf 70 (ouf of a maximum of 159) associated hifs in fhe 
TPC, a ;^^/ndf < 2 of fhe frack momenfum fif in fhe TPC and af leasf one associated hif in one of fhe fwo 
SPD layers. Wifh fhese frack selection criferia, fhe accepfance in rapidify for D mesons drops steeply 
fo zero for |yl+ 0.5 af low px and for |y|+0.8 af px;c;5 GeV/c. A px-dependenf fiducial accepfance cut 
was therefore applied on the D-meson rapidity, |y| < yfid(pT)> with yfid(pT) increasing from 0.5 to 0.8 in 
0 < px < 5 GeV/c according to a second order polynomial function and taking a constant value of 0.8 
for Px > 5 GeV/c. 

D+ candidates were filtered by applying kinematical cuts and geometrical selections on the decay 
topology, together with particle identification criteria. The selection criteria were tuned in each px 
interval and centrality class to have a good statistical significance of fhe signal, while keeping fhe 


5 









Ds production in Pb-Pb collisions 


ALICE Collaboration 


selection efficiency as high as possible. It was also checked that background fluctuations were not 
causing a distortion in the signal line shape by verifying that the D+-meson mass and its resolution 
were in agreement with the Particle Data Group (PDG) world-average value (1.969 GeV/c^ |[52]l l and 
the Monte Carlo simulation results, respectively. The resulting selection criteria depend on the transverse 
momentum of the candidate and provide a selection efficiency that increases with increasing pj. 

The main variables used to select the D+ decay topology were the decay length (L), defined as the 
distance between the primary and secondary vertices, and the cosine of the pointing angle (cos 6point)> 
which is the angle between the reconstructed D+ momentum and the line connecting the primary and 
secondary vertices. Additional selections were applied on the projections of decay length and cosine of 
pointing angle in the transverse plane xy (Lxy, cos in order to exploit the better resolution on the 

track parameters in that plane. A further cut was applied on L^y divided by its uncertainty {L^y /The 
three tracks were also required to have a small distance to the reconstructed decay vertex, by defining the 
variable dvenex as the square root of the sum in quadrature of the distances of each track to the secondary 
vertex. To further suppress the combinatorial background, the angles d*{n), i.e. the angle between the 
pion in the KK;r rest frame and the KK;r flight line in the laboratory frame, and P'(K), i.e. the angle 
between one of the kaons and the pion in the KK rest frame, were exploited. The cut values used for D+ 
mesons with 4 < pT < 6 GeV/c in the 0-10% centrality class were: L, Lxy > 500 /im, Lxy/Oi^^. > 7.5, 
cosOpoint > 0.94, cosOj^jj^j > 0.94, avertex < 400 pm, cos6*(7r) > 0.05 and |cos^ 0^(K)| < 0.9. Looser 
selection criteria were used for D+ selection at higher pj and in more peripheral events, due to the lower 
combinatorial background. 

In addition, to select D+ mesons decaying in the considered 071+ mode, with 0 —K^K+, candidates 
were rejected if none of the two pairs of opposite-charged tracks had an invariant mass compatible with 
the PDG world average for the 0 mass (1.0195 GeV/c^ ISH). The difference between the reconstructed 
K+K^ invariant mass and world-average 0 mass was required to be less than 4 MeV/c^ (a selection that 
preserves about 70% of the signal) for D+ candidates in the three px intervals considered in the 0-10% 
centrality class, while looser selections were used for semi-central events. 

Particle identification was used to obtain a further reduction of the background. Compatibility cuts were 
applied to the difference between the measured signals and those expected for a pion or a kaon. A track 
was considered compatible with the kaon or pion hypothesis if both its dE/Ax and time-of-flight were 
within 3 a from the expected values. Tracks without a TOE signal (mostly at low momentum) were 
identified using only the TPC information and requiring a 2a compatibility with the expected AE/Ax. 
Triplets of selected tracks were required to have two tracks compatible with the kaon hypothesis and 
one with the pion hypothesis. In addition, since the decay particle with opposite charge sign has to be a 
kaon, a triplet was rejected if the opposite-sign track was not compatible with the kaon hypothesis. This 
particle identification strategy preserves about 85% of the D+ signal. 

Eor each candidate, two values of invariant mass can be computed, corresponding to the two possible 
assignments of the kaon and pion mass to the two same-sign tracks. Signal candidates with wrong mass 
assignment to the same-sign tracks would give rise to a contribution to the invariant-mass distributions 
that could potentially introduce a bias in the measured raw yield of D+ mesons. It was verified, both in 
data and in simulations, that this contribution is reduced to a negligible level by the particle identification 
selection and by the requirement that the invariant mass of the two tracks identified as kaons is compatible 
with the 0 mass. 

The invariant-mass distributions of the D+ candidates (sum of D+ and D^^ candidates) are shown in Eig.[T] 
in the three px intervals for the 10% most central Pb-Pb collisions. The raw signal yields were extracted 
by fitting the invariant-mass distributions with a function that consists of the sum of a Gaussian term to 
describe the signal peak and an exponential function to describe the background. The fit was performed 
in the invariant-mass range 1.88 < M(KK7r) <2.1 GeV/c^ in all px intervals. The lower limit of 1.88 
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Centrality class 

Pt interval 
(GeV/c) 

raw 

S/B (3a) 


;i^^/ndf 

S/VS+ B (3a) 


4-6 

438±144 

0.02 

3.0 

27.4/18 

0-10% 

6-8 

117± 38 

0.10 

3.2 

17.5/18 


8-12 

89± 21 

0.38 

5.0 

26.5/18 

20-50% 

6-8 

197± 61 

0.07 

3.5 

9.9/21 


8-12 

52± 20 

0.29 

3.4 

17.9/21 


Table 2: Measured raw yields signal over background (S/B), statistical significance ( S/-\/S + B ) and 

;i;^/ndf of the invariant-mass fit for D+ and their antiparticles in the considered pj intervals for the 0-10% and 
20-50% centrality classes. The uncertainty on the raw yield is the statistical uncertainty obtained from the fit. 

GeV/c^ was chosen to exclude the contribution of D+ —)■ K^K+;r+ decays, BR = (0.265 ^qqq 9)% HH, 
which could give rise to a bump in the background shape for invariant-mass values around the D+ mass 
(1.870 GeV/c^) |[5^ . The mean values of the Gaussian functions in all the pj intervals are compatible 
within two times their uncertainty with the PDG world average for the D+ mass and the Gaussian widths 
are in agreement with the expected values from Monte Carlo simulations. 



Figure 1: Invariant-mass distributions of D+ candidates and charge conjugates in the three considered pj intervals 
in the 10% most central Pb-Pb collisions. 


In Table[2]the extracted raw yields of D+ mesons (sum of particle and antiparticle), defined as the integral 
of the Gaussian functions, are listed for the different pt intervals in both the considered centrality classes, 
together with the signal-over-background (S/B) ratios and the statistical significance (S/\/S + B). The 
background was evaluated by integrating the background fit functions in ±3a around the centroid of the 
Gaussian. 


4 Corrections 


The raw yields extracted from the fits to the invariant-mass distributions of D+ and candidates were 
corrected to obtain the production yields of prompt (i.e. not coming from weak decays of B mesons) D+ 
mesons. The pj-differential yield of prompt D+ was computed as 


AN 


D+ 


dpT 


/prompt (PT ) 


raw 


(pt) 


l.vKffld 


b|<0.5 


ApT BR'Aevt 2yfi(j(p'p) (AcC X £)prompt(PT) 


( 1 ) 


where ”'^(pt) are the values of the raw yields (sum of particles and antiparticles) reported in Table|2j 
which were corrected for the B-meson decay feed-down contribution (i.e. multiplied by the prompt 
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Figure 2: Left; Acceptance-times-efficiency for D+ mesons in the 10% most central Pb-Pb collisions. The 
efficiencies for prompt (solid lines) and feed-down (dotted lines) D+ mesons are shown. Also displayed, for 
comparison, the efficiency for prompt D+ mesons without PID selections (dashed lines). Right: Relative variation 
of the prompt D+-meson yield in the 0-10% centrality class as a function of the hypothesis on 
for the B feed-down subtraction approach based on Eq. ©. 

fraction /prompt)^ divided by the acceptance-times-efficiency for prompt D+ mesons, (Acc x £)prompt> and 
divided by a factor of two to obtain the charge (particle and antiparticle) averaged yields. The corrected 
yields were divided by the decay channel branching ratio (BR), the px interval width (Apx). the rapidity 
coverage (2yfid) and the number of analysed events (Aevt). 

The correction for the acceptance and the efficiency was defermined using Monfe Carlo simulafions. 
Pb-Pb collisions af = 2.76 TeV were simulated using the HIJING vL383 event generator |[54]| . 
Prompt and feed-down D+ (and signals were added with the PYTHIA v6.4.21 generator ll55l . In 
order to minimize the bias on the detector occupancy, the number of D mesons injected into each HIJING 
event was adjusted according to the Pb-Pb collision centrality. The px distribution of the generated D+ 
mesons in the 0-10% centrality class was weighted in order to match the shape measured for D° mesons 
in central Pb-Pb collisions ll42]| . For the 20-50% centrality class, the generated px distribution was 
defined based on FONLL perfurbafive QCD calculafions II561I57II mulfiplied by fhe nuclear modificalion 
faclor predicfed by fhe DAMPS parfonic fransporf model |[58]l . which reproduces fhe measured non- 
sfrange D-meson Raa in semi-cenfral collisions wifhin uncerfainfies ifT^ . 

The generated particles were fransporfed fhrough fhe AFICE defector using fhe GEANT3 1591 particle 
fransporf package fogefher wifh a defailed description of fhe geomefry of fhe apparafus and of fhe 
defector response. The simulafion was funed to reproduce fhe position and widfh of fhe inferacfion 
vertex disfribufion, fhe number of active elecfronic channels and fhe accuracy of fhe defecfor calibrafion, 
and fheir fime evolution wifhin fhe Pb-Pb dafa faking period. 

The efficiencies were evaluated in cenfralify classes corresponding fo fhose used in fhe analysis of fhe 
dafa in ferms of charged-parficle mulfiplicify, hence of defector occupancy. In fhe leff-hand panel of 
Fig.|2j fhe (Acc x e) values for prompf and feed-down D+ mesons wifh rapidify |y| < yfiu are shown for 
fhe 0-10% cenfralify class. The same figure shows also fhe (Acc x e) values for fhe case wifhouf fhe PID 
selections, demonsfrafing fhaf fhis selection is abouf 85% efficienl for fhe signal. 

The magnifude of (Acc x e) increases wifh increasing px, from 0.4% in fhe lowesf px interval up fo 
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2% in 8 < pt < 12 GeV/c. The (Acc x e) values for D+ from beauty-hadron decays are larger than 
those for prompt D+ by a factor of approximately 2.5-3.5 depending on pj, because the decay vertices 
of the feed-down D+ mesons are more displaced from the primary vertex and they are, therefore, more 
efficiently selected by the analysis cuts. The efficiency of fhe selections used in the centrality interval 
20-50% is higher by a factor of about two with respect to that in the most central events, because the 
smaller combinatorial background in semi-peripheral collisions allowed the usage of looser selections 
on the D+ candidates. 

The ratio of prompt to inclusive contributions in the D+-meson raw yield, /prompt> was evaluated using a 
procedure similar to the one adopted for the pp measurement Il45]l . The contribution of feed-down from B 
decays in the raw yield depends on pj and on the applied geometrical selection criteria. The feed-down 
contribution was estimated using the beauty-hadron production cross section from FONLL perturbative 
QCD calculations for pp collisions at = 2.76 TeV scaled by the average nuclear overlap function 
(Taa) in each centrality class, the B—)-D-i-A decay kinematics from the EvtGen package |[60]l and the 
Monte Carlo efficiencies for feed-down D+ mesons. The resulting sample of feed-down D+ mesons is 
composed of two contributions: about 50% of the feed-down originates from B^-meson decays, while the 
remaining 50% comes from decays of non-strange B mesons (B° and B+). A hypothesis on the nuclear 
modification factor of feed-down D+ mesons, was introduced to account for the different 

modification of beauty and chai'm production in Pb-Pb collisions and for the possible enhancement of 
the Bg over non-strange B-meson yield due to the effect of hadronisation via recombination 1^ . The 
fraction of prompt D+ yield was therefore computed in each pj interval as 


./prompt — 1 


N- 


feed-down raw 




= ( 1 ^) 


FONLL 


feed-down 


R 


ifeed-down (AcC X £)feed-down ' 2yfid Apx ‘ EE ‘ Aevt 


( 2 ) 


AA 


raw 


/2 


where (Acc x £)feed-down is the acceptance-times-efficiency for feed-down D+ mesons. To determine the 
central value of /prompt> it was assumed that the nuclear modification factors of feed-down and prompt 
D+ mesons were equal The resulting feed-down contribution is about 20-25% 

depending on the pj interval. To determine the systematic uncertainty the hypothesis was varied in the 
range 1/3 < < 3, as discussed in detail in Section[5l It should be noted that the central 

value and the range of the hypothesis on differ from those used for non-strange D 

mesons in Refs. lEiiiiiiia, owing to the unknown role of recombination in the beauty sector, which 
could enhance the ratio of B° over non-strange B mesons, and to the large fraction of feed-down D+ 
mesons originating from non-strange B-meson decays. 

The nuclear modification factor of D+ mesons was computed as 


Eaa(pt) 


dA^AA/dPT 
{TAA)<i(y^p /dpT 


(3) 


The values of the average nuclear overlap function, (Taa), for the considered centrality classes are 
reported in Table[I] The pT-differential cross section of prompt D+ mesons with |y| < 0.5 in pp collisions 
at = 2.76 TeV, used as reference for Raa, was obtained by scaling in energy the measurement at 
y/s = 1 TeV ll45l . The ratio of the cross sections from FONEE pQCD calculations ll57l at y/s = 2.76 
and 7 TeV was used as the scaling factor. Since FONEE does not have a specific prediction for D+ 
mesons, the cross sections of the D-meson admixture (70% of D** and 30% of D+) were used for the 
scaling. The theoretical uncertainty on the scaling factor was evaluated by considering the envelope of 
the results obtained by varying independently the factorisation and renormalisation scales and the charm 
quark mass, as explained in detail in Ref. |[62l . For D+ and D*+ mesons, the result of the scaling was 
validated by comparison with data Il63l . 
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0-10% centrality 

20-50% centrality 


Pj interval (GeV /c) 

Pj interval (GeV /c) 


4-6 

6-8 

8-12 

6-8 

8-12 

Raw yield extraction 

8% 

8% 

8% 

8% 

8% 

Tracking efficiency 

15% 

15% 

15% 

15% 

15% 

Selection efficiency 

20% 

20% 

20% 

20% 

20% 

PID efficiency 

7% 

7% 

7% 

7% 

7% 

MC pt shape 

2% 

1% 

1% 

1% 

1% 

Feed-down from B 






FONFF feed-down corr. 

-28 

110% 

-21% 

+ 

-20 

-25 

y^fed-down /^prompt 

+ 10% 
-22 

-> 


+ 11% 
-22/c 

+ 12% 

-24/c 

Centrality limits 

< 1% 

< 1% 

Branching ratio 

4.5% 


Table 3: Relative systematic uncertainties on px-differential yields of prompt D+ mesons in Pb-Pb collisions for 
the two considered centrality classes. 

5 Systematic uncertainties 

The systematic uncertainties on the prompt D+-meson yields in Pb-Pb collisions are summarised in 

Tables 

The systematic uncertainty on the raw yield extraction was estimated from the distribution of the results 
obtained by repeating the fit to the invariant-mass spectra varying i) the fit range and ii) the probability 
distribution functions used to model the signal and background contributions. In particular, a second 
order polynomial function was used as an alternative functional form to describe the background. The 
signal line shape was varied by using Gaussian functions with mean and width fixed to the world-average 
D+ mass and to the values expected from Monte Carlo simulations, respectively. Furthermore, the raw 
yield was also extracted by counting the entries in the invariant-mass distributions after subtraction 
of the background estimated from a fit to the side bands of the D+ peak. In case of fitting in an 
extended mass range, it was verified that the effect on the D+ yield due to the possible bump produced 
in the candidate invariant-mass distribution by D+ —K^K+;r+ decays was negligible. An 
additional test was performed by fitting the D+ candidate invariant-mass distribution after subtracting 
the background estimated by coupling a pion track with K+K^ pairs having an invariant mass in the side 
bands of the (j) peak. The uncertainty was estimated to be 8% in all pj intervals. 

The contribution to the measured yield from D+ decaying into the K^K+;r+ final state via other resonant 
channels (i.e. not via a (p meson) was found to be negligible, due to the much lower selection efficiency, 
as discussed in Ref. na. 

Ofher contributions to the systematic uncertainty originate from the imperfect implementation of the 
detector description in the Monte Carlo simulations, which could affect the particle reconstruction, the 
D+ selection efficiency, and the kaon and pion identification. 

The systematic uncertainty on the tracking efficiency (including the effect of the track selection) was 
estimated by comparing the efficiency (i) of track finding in the TPC and (ii) of track prolongation from 
the TPC to the ITS between data and simulations, and (iii) by varying the track quality selections. The 
estimated uncertainty is 5% per track, which results in 15% for the three-body decay of D+ mesons. 

The effect of residual discrepancies between data and simulations on the variables used to select the D+ 
candidates was estimated by repeating the analysis with different geometrical selections on the decay 
topology and varying the cut on the compatibility between the K+K^ invariant mass and the (j) mass. A 
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systematic uncertainty of 20% was estimated from the spread of the resulting corrected yields. 


The systematic uncertainty induced by a different efficiency for particle identification in data and sim¬ 
ulations was estimated by comparing the corrected D+ yields obtained using different PID approaches, 
testing both looser and tighter cuts with respect to the baseline selection described in Section |4l Due to 
the limited statistical significance, an analysis wifhouf PID selection could nof be carried ouf. Such a fesf 
was performed in fhe analysis of D*’ (—)■ K^;r+), D+ (—;■ K^;r+7r+) and D*+ (—;■ D‘^7r+) and a 5% uncer- 
fainfy was estimated for fhe case of 3 a cufs on dE/dx and fime-of-llighl signals, which correspond fo 
fhe loosesl selecfions fhaf could be tested for fhe D+. Based on all fhese checks a sysfemafic uncertainly 
of 7% on fhe PID selecfion efficiency was eslimaled. 

The efficiency is also sensifive fo differences belween fhe real and simulated D+ momenfum disfribufions. 
The effecf depends on fhe widfh of fhe pj intervals and on fhe variation of fhe efficiency wilhin Ihem. A 
sysfemafic uncertainly was defined from fhe relalive difference among fhe efficiencies oblained using 
differenl px shapes for fhe generated D+ mesons, namely fhe measured dA/dpx of D*^ mesons in 
cenlral Pb-Pb collisions, fhe px shape predicfed by FONLL pQCD calculafions wilh and wifhouf fhe 
nuclear modification predicfed by fhe DAMPS parfonic fransporf model. The resulting confribufion fo 
fhe systematic uncertainly was found fo be 2% for fhe momenfum inlerval 4 < px < 6 GeV/c, where fhe 
selecfion efficiency is slrongly px dependenl, and 1% al higher px. 


The sysfemafic uncerlainfy due fo fhe subfraclion of D+ mesons from B-meson decays was eslimaled 
following fhe procedure described in Ref. ifTTI . The confribufion of fhe uncerfainfies inherenf in fhe 
FONLL perlurbalive calculalion was included by varying fhe heavy-quark masses and fhe factorisation 
and renormalisalion scales, pp and Pr, independenlly in fhe ranges 0.5 < pp/mx < 2, 0.5 < Pr/wix < 2, 
wilh fhe consfrainl 0.5 < pp/pR < 2, where nix = ■\JPj + wiq- Furlhermore, fhe prompf fracfion oblained 
in each px interval was compared wilh fhe resulls of a differenl procedure in which fhe FONLL cross 
secfions for prompf and feed-down D mesons and Iheir respeclive Monte Carlo efficiencies were fhe 
inpuf for evalualing fhe correction facfor 


V' 

1 I (AcC X g)feed-down \ dvdpx ) feed-down ^AA 

(AcC X £)pro„,pt ■ FONLL • 

\ V / prompt / 


(4) 


Since FONLL does nof have a specific prediction for D+ mesons, four differenl approaches were used 
to compute fhe predicted px shapes of promplly produced D+, (d^a/dydpx)pj*^J)|pj'", as explained in 
delail in Ref. Il45]| : (i) FONLL prediction for fhe admixlure of charm hadrons; (ii) FONLL prediclion 
for D*+ mesons (fhe D*+ mass being close fo fhaf of fhe D+); (iii) FONLL prediclion for c quarks and 
fragmenlalion funclions from ll64ll wilh parameler r = (wid — n^c )/wid and being fhe masses of fhe 
considered D-meson species and of fhe c quark, respeclively); (iv) FONLL prediction for c quarks and 
fragmenlalion functions from llMl wilh parameter r = 0.1 (as used in FONLL calculafions) for all meson 
species. In fhe latter fwo cases, fhe D*+ mesons produced in fhe c quark fragmenlalion were made to 
decay wilh PYTHIA and fhe resulting D+ were summed fo fhe primary ones fo obfain fhe prompf yield. 
The sysfemafic uncerlainfy due fo fhe B feed-down subfraclion was finally evaluated as fhe envelope of 
fhe resulls oblained wilh fhe fwo melhods, namely Eq. Q and (@1), when varying fhe FONLL paramelers 
and fhe c—;■ D+ fragmenlalion function used fo determine (d^a/dydpx)p^*^^p|'" in Eq. dUl. 

The confribufion due fo fhe differenl nuclear modification faclor of prompf and feed-down D+ mesons 
was eslimaled by varying fhe hypolhesis on “ the range 1/3 < < 3 

for bolh feed-down subfraclion melhods. The varialion of fhe hypolhesis is molivafed by fhe combined 
effecf on fhe Raa of (i) the differenl energy loss of charm and beaufy quarks in fhe QGP, as predicfed 
by energy loss models and supporled by experimenfal dala on D meson and non-prompf J/t/r Raa af fhe 
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px interval (GeV/ c) 


4-6 

6-8 

8-12 

Data systematics in pp 

26% 

25% 

29% 

Eeed-down from B 

+ 4% 
-17 


^17% 

-yi-scaling of the pp reference 

1 ' 4 % 


+ «% 

Normalisation 


3.5% 


Branching ratio 


4.5% 



Table 4: Relative systematic uncertainties on the pp reference cross section. The row labeled ‘Data systematics’ 
reports the sum in quadrature of the contributions due to raw yield extraction, tracking efficiency, selection 
efficiency, PID efficiency, MC px shape and ’other resonant channels’ from Ref. Il45l . 

LHC II111I131I421I651I66II ; (ii) the possibly different contribution of coalescence in charm and beauty quark 
hadronisation, leading to a different abundance of D+ and B° mesons relative to non-strange mesons; and 
(iii) the possibly different modulation of D and B spectra due to radial flow. The resulting uncertainty 
for the case of B feed-down subtraction approach based on Eq. Q is shown in the right-hand panel of 
Eig.|2]for the thr'ee pi intervals in the 0-10% centrality class. 

The Pb-Pb data are also affected by a systematic uncertainty on the determination of the limits of the 
centrality classes, due to the 1.1% relative uncertainty on the fraction of the total hadronic cross section 
used in the Glauber fit lISTTl . This contribution was estimated from the variation of the D-meson dA/dpx 
when the limits of the centrality classes are shifted by ±1.1%. The resulting uncertainty, which is 
common to all pj bins, is less than 1% for both the 0-10% and the 20-50% centrality classes. 

Einally, the 4.5% uncertainty on the branching ratio ll52l was considered. 

In the calculation of the /?aa> the uncertainties on the reference cross section for pp collisions, the Pb-Pb 
yields, and the average nuclear overlap function were considered. 

Eor the pp reference, the uncertainties on the measurement at ^/s = 7 TeV, described in Ref. ll45l and 
those due to the EONLL-based scaling to ^/s = 2.76 TeV, described in Section |4l were summed in 
quadrature. The contributions to the systematic uncertainty on the pp reference cross section are reported 
in Table |4] 

The uncertainties on the pp reference were added in quadrature to those on the Pb-Pb prompt D+ yields, 
described above, except for the BR that cancels out in the ratio and the feed-down contribution deriving 
from EONLL uncertainties, that partly cancels in the ratio. This contribution was evaluated by comparing 
the /?aa values obtained with the two methods for feed-down correction of Eq. (|2l) and (lUl and with the 
different heavy quark masses, fragmentation functions, factorisation and renormalisation scales used in 
EONLL. In this study, these variations were done simultaneously for the Pb-Pb yield and for the pp 
reference cross section, so as to take into account the correlations of these sources in the numerator and 
denominator of /?aa- 

Einally, the /?aa normalisation uncertainty was computed as the quadratic sum of the 3.5% pp normal¬ 
isation uncertainty Il45]l . the contribution due to the 1.1% uncertainty on the fraction of hadronic cross 
section used in the Glauber fit discussed above, and the uncertainty on (Taa), which is of 3.2% and 3.7% 
for the 0-10% and 20-50% centrality classes, respectively. 

6 Results 

The transverse momentum distributions dV/dpx of prompt D+ mesons in Pb-Pb collisions are shown in 
Eig-El for the 0-10% and 20-50% centrality classes. The yields reported in Pig.|3]refer to particles only. 
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Figure 3: Transverse momentum distributions dN/dpj of prompt D+ mesons in the 0-10% (left panel) and 
20-50% (right panel) centrality classes in Pb-Pb collisions at = 2.76 TeV. Statistical uncertainties (bars), 
systematic uncertainties from data analysis (empty boxes) and systematic uncertainties due to beauty feed-down 
subtraction (shaded boxes) are shown. The reference pp distributions (7 aa) dcr/dpx are shown as well. 


since they were computed as the average of particles and antiparticles under the assumption that the 
production cross section is the same for D+ and D^. The vertical error bars represent the statistical 
uncertainties. The symbols are positioned horizontally at the centre of each pj interval, with the 
horizontal bars representing the width of the pj interval. The systematic uncertainties from data analysis 
are shown as empty boxes around the data points, while those due to the B feed-down subtraction, 
which include the contributions of the FONLL uncertainties and of the variation of the hypothesis on 
are displayed as shaded boxes. The normalisation uncertainties are reported as text on 

the figures. 

The pT-differential yields measured in Pb-Pb collisions are compared to the reference yields in pp 
collisions at the same energy, scaled by the nuclear overlap function (T^a), reported in Table [B The 
pp reference at ^/s = 2.76 TeV is obtained by scaling the cross section measured at 7 TeV as described in 
Section m A clear suppression of the D+-meson yield in the 10% most central Pb-Pb collisions relative 
to the binary-scaled pp yields is observed in the highest pj interval (8 < px < 12 GeV/c). In the 20- 
50% centrality class, an indication of suppression is found in 8 < px < 12 GeV/c. At lower px, in both 
centrality classes, it is not possible to conclude on the presence of a suppression of the D+-meson yield 
in heavy-ion collisions with respect to the pp reference. 

The nuclear modification factor Raa of prompt D+ mesons was computed from the dA/dpx distributions. 
The results are shown as a function of px in the left-hand panel of Fig.|4]for the two centrality classes. The 
vertical bars represent the statistical uncertainties, the empty boxes are the total px-dependent systematic 
uncertainties described in Section |5j except for the normalisation uncertainty, which is displayed as a 
filled box at Raa = 1 ■ A suppression by a factor of about three of the D+-meson yield in Pb-Pb collisions 
relative to the binary-scaled pp cross section is observed in the highest px interval (8 < px < 12 GeV/c) 
for the 10% most central collisions. A smaller suppression (by a factor of about two) is measured 
in the 20-50% centrality class in 8 < px < 12 GeV/c, even though with the current uncertainties no 
conclusions can be drawn on the centrality dependence of the D+-meson nuclear modification factor 
at high Px. Since no significant modification of the D+-meson production relative to binary-scaled pp 
collisions is observed in p-Pb reactions in the px range considered here |[T4]| . the substantial suppression 
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Figure 4: Left: Raa of prompt D+ mesons in the 0-10% and 20-50% centrality classes as a function of pj. 
For the 20-50% case, the symbols are displaced horizontally for visibility. Right: Raa of prompt D+ mesons 
compared to non-strange D mesons (average of D+ and D*+ 1421) in the 0-10% centrality class. Statistical 
(bars), systematic (empty boxes), and normalisation (full box) uncertainties are shown. 
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of the D+-meson yield at high pj in Pb-Pb collisions cannot be explained in terms of initial state effects, 
but it is predominantly due to strong final-state effects induced by the hot and dense partonic medium 
created in the collisions of heavy nuclei. At lower pj the central values of the measurement show a larger 
Raa, however the large statistical and systematic uncertainties do not allow to draw a conclusion on the 
Pj dependence of the D+ nuclear modification factor. 

The Raa of prompt D+ mesons in the 10% most central collisions is compared in the right-hand panel 
of Fig. |4] to the average nuclear modification factor of D*’, D+ and D*+ mesons measured in the same 
centrality class ll42]| . This comparison is meant to address the expected effect of hadronisation via quark 
recombination in the partonic medium on the relative abundances of strange and non-strange D-meson 
species. In the three pj intervals, the values of the D+-meson Raa are higher than those of non-strange 
D mesons, although compatible within uncertainties. Even considering that a part of the systematic 
uncertainty is correlated between strange and non-strange D mesons, the current uncertainties do not 
allow a conclusive statement on the expected enhancement of the D+-meson yield relative to that of 
non-strange D mesons in heavy-ion collisions. 

An alternative approach to study the predicted modification of the chai'm-quark hadronisation in the 
presence of a QGP is to compare the ratios between the measured yields of D+ and D°(D+) mesons 
in Pb-Pb and pp collisions. This comparison is shown in Fig. [5] for the 10% most central Pb-Pb 
collisions. In the left-hand panel the D+/D° ratio is displayed, while the right-hand panel shows the 
ratio D+/D+. The ratios D+/D° and D+/D+ in pp collisions are taken from the measurements at 
y/s = 1 TeV lldsfl No strong dependence on the collision energy is expected (see ll45l and references 
therein). In the evaluation of the systematic uncertainties on the D-meson yield ratios, the sources of 
correlated and uncorrelated systematic effects were treated separately. In particular, the contributions of 
the yield extraction, topological selection efficiency and PID efficiency were considered as uncorrelafed 
and summed in quadrafure. The uncerfainfy on fhe fracking efficiency cancels complefely in fhe ratios 
befween producfion cross sections of meson species reconsfrucfed from fhree-body decay channels (D+ 

^ The values from Ref. m were re-computed with the most recent value for the branching ratio of the 
D+ —> ((i;r+ K^K+;r+ decay chain, which is 2.24% (52), while it was 2.28% at the time of the pp publication. 
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Figure 5: Ratios of prompt D-meson yields (D+ZD** and D+/D+) as a function of pj in the 10% most central 
Pb-Pb collisions at = 2.76 TeV compared to the results in pp collisions at ^ = 1 TeV. Statistical (bars) and 
systematic (boxes) uncertainties are shown. 
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and D+), while a 5% systematic uncertainty (4% in the pp case) was considered in the ratio to the D° 
yields, which are reconstructed from a two-particle final state. To propagate the uncertainty due to the B 
feed-down subtraction, the contribution of the FONLL cross section was treated as completely correlated 
among the D-meson species. It was estimated from the spread of the D-meson yield ratios obtained by 
varying the factorisation and renormalisation scales and the heavy-quark mass in FONLL coherently for 
the three meson species. The contribution due to the hypothesis on was considered as 

uncorrelated between D+ and non-strange D mesons and summed in quadrature. The difference between 
the D+/D** ratios in pp and in central Pb-Pb collisions is of about 1 a of the combined statistical and 
systematic uncertainties in both the two lowest pi intervals, 4 < pj <6 GeV/c and 6 < pj <S GeV/c. 
An enhancement of Dg/D ratios in heavy-ion collisions is predicted if recombination contributes to charm 
quark hadronisation in the QGP However, considering the current level of experimental uncertainties, 
no conclusion on charm-quark hadronisation can be drawn from this first measurement of D+-meson 
production in Pb-Pb collisions. 

In the framework of the Statistical Hadronisation Model II391I671I68II . the py-integrated ratios of D-meson 
abundances for a chemical freeze-out temperature T = 156 MeV (as extracted from fits to the measured 
abundances of light-flavour hadrons |[69l ) and vanishing baryo-chemical pofenfial, are expecfed fo be 
D+/D° = 0.338 and D+/D+ = 0.830, which are higher by a factor of abouf fwo wifh respecf fo fhe 
values calculafed for pp collisions af LHC energies ll43l . 

In Fig. fhe measured Raa of non-sfrange D mesons and of D+ are compared fo fhe predicfion of 
fhe TAMU model EtIIMII . Among fhe several models available for open charm producfion in heavy- 
ion collisions, TAMU is fhe only one providing a quanfifafive predicfion for fhe D+-meson nuclear 
modificalion faclor. This is a heavy-quark fransporf model based on heavy-quark diffusion, implemenfed 
via simulafions based on fhe relafivisfic Langevin equafion, in a hydrodynamically expanding medium. 
The interacfions of fhe charm quarks wifh fhe medium are modeled including only elasfic processes, 
which are assumed fo govern fhe heavy-quark scattering amplifudes af low and infermediafe momenfa. 
The heavy-quark fransporf coefficienfs are calculafed wifhin a non-perfurbafive T-mafrix approach, 
where fhe interacfions proceed via resonance formafion fhaf fransfers momenfum from fhe heavy quarks 
fo fhe medium consfifuenfs. The hadronisafion of charm quarks is performed via recombinafion wifh 
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Figure 6: Raa of prompt D+ and non-strange D mesons (average of D**, D+ and D*+) in the 0-10% centrality 
class compared to predictions of the TAMU model ED- The bands shown for the TAMU predictions encompass 
the charm-shadowing uncertainty. 


thermalized up, down and strange quarks. The remaining charm quarks are converted to hadrons using 
the vacuum fragmentation functions from Il64ll and fragmentation fractions /(c —)■ D) from PYTHIA. 
This model predicts an enhancement of the D+ over the non-strange D-meson Raa at low px as a 
consequence of the recombination of charm quarks with thermally equilibrated strange quarks in the 
QGP. At higher pj, where the dominant hadronisation mechanism is fragmentation, similar Raa values 
are predicted for the different D-meson species. The model describes the measured D+-meson nuclear 
modification factor within uncertainties and at low pj provides also a reasonable description of non- 
strange D-meson Raa- The measured suppression of non-strange D mesons is underestimated at higher 
Pi, where the contribution of inelastic processes (gluon radiation), which are missing in this transport 
calculation, is expected to play a major role. 

7 Summary 

The production of D+ mesons was measured for the first time in heavy-ion collisions. The measurement 
was carried out on a sample of Pb-Pb collisions at = 2.76 TeV in two centrality classes, namely 
0-10% and 20-50%. 

The results for the 10% most central collisions indicate a substantial suppression {Raa ~ 0.3) of the 
production of D+ mesons at high pj (8 < px < 12 GeV/c) with respect to the expectation based on the 
pp cross section scaled by the average nuclear overlap function. The observed suppression is compatible 
with that of non-strange D mesons and can be described by models including strong coupling of the 
charm quarks with the deconfined medium formed in the collision. 

At lower momenta (4 < px < 8 GeV/c), the values of the D+-meson nuclear modification factor are 
larger than those of non-strange D mesons, although compatible within uncertainties. This result provides 
a possible hint for an enhancement of Dj/D ratio, which is expected if the recombination process 
significantly contributes to the charm quark hadronisation in the QGP. 

The precision of the measurements will be improved using the larger data samples of Pb-Pb collisions 
that will be collected during the ongoing LHC Run-2. The larger sample size will allow us to observe 
the D+ signal with less stringent selections, thus reducing the systematic uncertainty on the efficiency 
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correction. In addition, the higher Pb-Pb collision centre-of-mass energy will reduce the impact of the 
-y/J-scaling of the pp reference. This will open the possibility to exploit the measurement of D+-meson 
production in heavy-ion collisions to assess the recombination effects in the charm-quark hadronisation 
and to provide further constraints to models describing the coupling of heavy quarks with the medium. 
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